Ab initio study of structural stability of small 3d late transition metal clusters 
Interplay of magnetization and hybridization 
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Using first-principles density functional theory based calculations, we analyze the structural sta- 
bility of small clusters of 3d late transition metals. We consider the relative stability of the two 
structures - layer-like structure with hexagonal closed packed stacking and more compact structure 
of icosahedral symmetry. We find that the Co clusters show an unusual stability in hexagonal sym- 
metry compared to the small clusters of other members which are found to stabilize in icosahedral 
symmetry based structure. Our study reveals that this is driven by the interplay between the mag- 
netic energy gain and the gain in covalency through s-d hybridization effect. Although we have 
focused our study primarily on clusters of size 19 atoms, we find that this behavior to be general 
for clusters having sizes between 15 and 20. 
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I. INTRODUCTION 

Atomic clusters of nanometer size have attracted a spe- 
cial attention of the present day research due to their 
applications in the field of optoelectronic catalysis;^ 
data-storage?^ sensors^"— etc. The first step to the theo- 
retical study of the properties of clusters is the determina- 
tion of the minimum energy structures. The equilibrium 
minimum energy structures of small clusters often prefer 
compact geometriesi^ like icosahedral or cub-octahedral 
symmetry based structures. It has also been found in sev- 
eral cases that the deformed three dimensional sections of 
the face centered cubic (fee) or hexagonal closed packed 
(hep) lattice appear as degenerate energy state or closely 
lying isomersJ^ However, depending upon the local sym- 
metry, they exhibit profoundly different properties. For 
example, planner gold clusters exhibit outstanding cat- 
alytic activity compared to their bulk counterpart of fee 
symmetry^i^ bi-layer Ru-nanoclusters exhibit significant 
chemical activity towards H2O splitting compared to Ru- 
clusters of hep symmetryJ^ Similarly, the dependence of 
magnetic behavior of the Pd clusters on cluster symmetry 
is found to be significant All of these indicate that the 
determination of the local symmetry is an unavoidable 
part in a cluster calculation. 

In this article, we have performed a first principles 
based analysis to understand the structural trend of tran- 
sition metal clusters. Transition metal clusters demand 
special attention because of their fascinating magnetic 
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the dependence of the equilibrium struc- 



ture on magnetisroii as well as their potential biomedical 
applications ji^"— We focus our attention only to the 3d 
late transition metal clusters. Among the 3d late tran- 
sition metal elements, Mn has a half- filled d-level, while 
the others have more than half-filled d-level. Consider- 



ing the earlier studies on structure of the 3d late transi- 
tion metal clusters, it is seen that the small Mn and Fe 
clusters generally prefer a compact icosahedral growth 
pattern as has been shown by the first principles cal- 
culations for the Mn clusters^i"— and Fe clusters 
For relatively less magnetic Ni clusters and nonmag- 
netic Cu-clusters, the first principles calculations and also 
some experimental evidences indicate mainly icosahedral 
growth pattern<^"— 1^"— Some recent calculations^"— 
also highlight non-icosahedral or amorphous structural 
pattern for small clusters of coinage metals like Cu, Ag 
and Au. On the other hand, the ferromagnetic small 
Cobalt cluster is quite different from the other members 
of the 3d late transition metal series, particularly Mn 
and Fe clusters. The small Co„ clusters rather prefer rel- 
atively non-compact layer-like structures. In our recent 
work^ using first-principles density functional study, we 
have shown a clear hexagonal growth pattern for the 
small Co„ clusters (15<n<20). Hexagonal symmetry 
based structures in this size range consist of 3-planes 
with hep stacking. Also recently, this layer-like structures 
of the small Co„ clusters (13<n<23) has been reported 
by Gong et a^t^ using density functional calculation. 
Experimental work on small Co„ clusters (n<50))^"— 
though is unable to give any definitive conclusion, indi- 
cates non-icosahedral packing too. 

It is therefore curious why the small cobalt clusters 
prefer a hep growth pattern with layer-like stacking, 
while the clusters of the other 3d late transition metal 
elements apparently prefer a more compact icosahedral 
growth pattern. In order to have an understanding on 
this issue, we have chosen these two close packed struc- 
tures with hep and icosahedral symmetries as the starting 
guesses and allowed them to relax under the assumption 
of collinear magnetic ordering. We have studied the rel- 
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ative stability between these two symmetry based struc- 
tures in terms of energetics, structural and electronic 
properties. We have carried out our study for the entire 
series of 3c? late transition metal clusters i.e Mn„, Fe„, 
Co„, Ni„ and Cu„. Our study reveals that the contrast- 
ing behavior of the stability of Co-clusters compared to 
the other members, arises due to the interesting interplay 
of the effects of magnetization and hybridization. 



II. COMPUTATIONAL DETAILS 

We employed density functional theory within plane 
wave pseudo potential method as implemented in Vienna 
ab initio simulation package<^ We used the projected 
augmented wave pseudo-potentials^i^ and the Perdew- 
Bruke-Ernzerhof exchange-correlation functional^ of the 
generalized gradient approximation (GGA). The pseudo 
potentials for the transition metal elements studied in 
this work, were generated considering the 3o? and 4s elec- 
trons as the valence electrons. The energy cut-off was 335 
eV for the cluster calculation of each transition metal. 
We did both spin-polarized and non-spin polarized cal- 
culations at the r point of the Brillouin zone. Geometry 
optimizations have been performed using the conjugate 
gradient and the quasi-Newtonian methods until all the 
force components were less than a threshold value of 0.01 
eV/A. For the cluster calculation, a simple cubic super- 
cell was used with periodic boundary conditions, where 
two neighboring clusters were kept separated by around 
12 A vacuum space, which essentially makes the interac- 
tion between cluster images negligible. To determine the 
magnetic moment of the minimum energy structure in 
spin-polarized calculations, we have explicitly considered 
all the possible spin multiplicities for each structure un- 
der the approximation of collinear atomic spin arrange- 
ments. 



III. RESULTS AND DISCUSSIONS 

To make our conclusion regarding structural and elec- 
tronic properties of the clusters more general, we have 
chosen the cluster size of 19 atoms for each 3d late tran- 
sition metal cluster considered in this study, instead of 
cluster size of 13-atoms, which is the first geometric 
magic size for the icosahedral symmetry based structure 
of the most transition metal clusters. Fig. [T] shows the 
cluster structures of 19 atoms with hep and icosahedral 
symmetries. A 19-atoms cluster structure with hep sym- 
metry, can be viewed as stacking of three planes contain- 
ing 6, 7 and 6 atoms respectively in each of these planes. 
On the other hand, a 19-atoms cluster structure with 
icosahedral symmetry can be thought of two interpene- 
trating 13-atoms icosahedrons. As seen from FiglU a 19- 
atoms cluster structure of hep symmetry appears more 




FIG. 1: (Color online) Cluster structures of 19 atoms with 
(a) hexagonal closed packed and (b) icosahedral symmetries. 
These two competing structural symmetries have been consid- 
ered in this work to determine the minimum energy structure 
for each Xig cluster, X=Mn, Fe, Co, Ni and Cu. 



open in the sense of lesser value of average co-ordination 
of atoms as well as more layer-like structure compared to 
the 19-atoms icosahedral structure. In the first step of 
our optimization procedure, we have started with ideal 
hexagonal and icosahedral structure for 19 atom clusters 
and have optimized them. In the second step, we have 
randomly displaced few atoms in the optimized struc- 
ture, obtained in the first step, and have reoptimized to 
get the "final optimized" structure. The second step has 
been carried out considering all possible collinear spin ar- 
rangements of the atomic spins in each Xig (X = Mn, Fe, 
Co, Ni, Cu) cluster within the spin-polarized calculation. 
The optimized clusters do not have perfect hexagonal or 
icosahedral symmetry but are heavily distorted. 

First, to analyze the optimized structures, we de- 
fine the average nearest-neighbor bond length as (r) = 
^ J2i>j ^ijj where is the bond distance between the 
j-th and i-th atoms, and nt is the number of such bonds. 
In cluster calculation, we considered that the two atoms 
are said to be bonded if their inter atomic distance is 
within 2.75 A, which is larger than any of the nearest 
neighbor bulk bond-lengths of these 3c? late transition 
metal elements. Fig. [2] shows the plot of the average 
nearest-neighbor bond-lengths of the optimized struc- 
tures of Xig clusters for both the symmetries. It is seen 
that the average nearest-neighbor bond-lengths for the 
hep structures are consistently smaller than those of the 
icosahedral counterpart in agreement with the previous 
studyr^ which indicates that the net attraction of nu- 
cleus on outer shell electrons is effectively more for the 
hep symmetry based structure. As the c?-shell gets filled 
one by one electron from Mn Fc — ^ Co — s- Ni ^ Cu, 
the ion-electron interaction gets stronger, which increases 
the binding. On the other hand, electron-electron repul- 
sion also increases, which starts to downplay the gain 
in electron-ion attraction. On top of this effect, the in- 
creased atom-centered magnetic moments also play a sig- 



TABLE I: Total binding energies and total magnetic moments 
of the minimum energy structure of hep and icosahedral sym- 
metries for each X19 cluster, X—Mn, Fe, Co, Ni and Cu in 
spin-polarized calculations. 



Clusters 


Binding energy (eV) 


Magnetic moment (/is) 




hep 


icosa 


hep 


icosa 


Mni9 


43.87 


45.12 


15 


19 


Fei9 


64.35 


66.26 


58 


58 


C019 


72.01 


70.80 


39 


37 


Nii9 


65.13 
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47.04 


46.95 









nificant role, specially for the niembers left of Co along 
C019 to Fei9 to Mnig. 

In Fig. 3, we further show the average co-ordination 
number plotted for the optimized structures of Xig clus- 
ters in both symmetries. We find that average co- 
ordination for hep symmetry based structure is less than 
that of icosahedral symmetry based structure, giving rise 
to a more open geometry although average bond length is 
smaller for hep based structure compared to icosahedral 
structure (cf. Fig. 2). 

Binding energy for each X19 cluster is calculated as 
Eb{Xiq) = [19E{X)-E{Xig)], where E(X) and E(Xi9) 
are the total energy of an isolated X atom and that of 
a X19 cluster respectively. In such a definition, a posi- 
tive sign in Eb corresponds to binding. Table I shows 




FIG. 2: (Color online) Plot of the average nearest-neighbor 
bond lengths (r) (see text) for the optimized hep and icosa- 
hedral 19-atoms clusters of Mn, Fe, Co, Ni and Cu. Blue 
(dark) squares correspond to the data points for the hep struc- 
ture and orange (light) circles for data points of the icosahe- 
dral structure of each X19 cluster in spin-polarized calcula- 
tion. The corresponding bulk-values have been shown with 
empty triangles. The atomic valance electronic configuration 
for each element has been marked at the top of the figure. 
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FIG. 3: (Color online) Plot of the average co-ordination num- 
ber for the optimized hep and icosahedral 19-atoms clusters 
of Mn, Fe, Co, Ni and Cu. 



the total binding energy and the total magnetic moment 
of the optimized hep and icosahedral structures of each 
X19 cluster. It is seen that the icosahedral symmetry 
based structure is stabler than the hep symmetry based 
structure for Mnig and Fcig clusters. Conversely, the hep 
symmetry based structure is energetically more favorable 
than the icosahedral structure for the Co 19 cluster, while 
both the structures are almost degenerate for the Niig 
and Cui9 clusters within the accuracy of our calculations. 
Analyzing the atomic spin orientations in the optimized 
structures of both the symmetries of each X19 cluster, 
we found that the Mn-Mn interactions within the Mnig 
cluster is mostly antiferomagnetic for both the optimal 
hep and optimal icosahedral phases as mentioned by the 
earlier worksi^i^^i On the other hand, each of the Fcig, 
Coig and Niig clusters is ferromagnetic for either of the 
two structural symmetries, with decreasing total mag- 
netic moment because the atom centered magnetic mo- 
ments decrease as one goes along Feig — > Coig — > Niig. 
Fig. m shows the distribution of atomic magnetic mo- 
ments of each of the Xig clusters derived from both the 
symmetries together with the corresponding optimized 
structures. Note that bulk Mn is also anti-ferromagnetic 
and bulk Fe, Co, Ni are ferromagnetic (with magnetic 
moment per atom 2.2 ^13 for Fe,i^i^ 1.72 ^13 for Co^ 
and 0.616 /is for Ni^). The Cuig cluster is found to be 
non- magnetic with zero magnetic moment. 

To have a visual representation, we have plotted the 
binding energy per atom for the optimal hep and op- 
timal icosahedral structures of each Xig cluster in Fig. 
El To understand the effect of magnetism on stability, 
we have also performed the non-spin polarized calcula- 
tion for each Xig cluster. The binding energies for the 
optimal hep and icosahedral structures of X19 clusters 
in non-spin polarized calculation are also shown in Fig. 
[5] (right panel) with shaded bars. Interestingly, non-spin 



4 



icosahedron 








FIG. 4: (Color online) Structures and atomic magnetic moment distribution in the optimized hep and optimized icosahedral 
structures of Mnig, Feig, C019 and Niig clusters in spin polarized calculation. For Mnig clusters, gray color represents up or 
positive and red (dark gray) color represents down or negative magnetic moment. For each of Feig, C019 and Niig clusters, the 
atoms are ferromagnetically coupled, each with positive or up magnetic moments and are, therefore, represented by same color 
(gray). The individual atomic magnetic moments of the constituent atoms in each optimized structure, have been represented 
by bar plots, length of the bars corresponding to the magnitude of atomic moments. The numbering of the atoms in each 
structure is indicated in the plots of the structures shown as insets. 



polarized calculation shows that the Mnig, Feig and C019 
clusters all stabilize in hep symmetry based structures, 
while both the structures are again degenerate for the 
Niig and Cuig clusters. From Fig. [51 it is also clearly seen 
that the magnetic phase always has the higher binding 
energy for both the structures of each Xig cluster, indi- 
cating that the magnetic phase is the stable phase for 
both the structures. Only in case of the Cuig cluster, 
the binding energy of each structure is same for both the 



spin polarized as well as non-spin polarized calculations, 
indicating that the Cuig cluster is nonmagnetic. From 
the binding energy plot, we therefore conclude that the 
magnetism switches the stable phase from hep to icosa- 
hedron in case of the Mnig and Fcig clusters, while the 
magnetism enhances further the stability of the hep phase 
for the Coig cluster. For the Niig and Cuig clusters, the 
effect of magnetism is small and both the hep as well as 
icosahedral symmetry based structures are almost degen- 




FIG. 5: (Color online) Plot of binding energy per atom of 
the minimum energy hep and icosahedral structures for each 
Xi9 {X— Mn, Fe, Co, Ni, Cu) cluster in both spin polarized 
calculation (left) and non-spin polarized calculation (right). 
The blue (dark) vertical bars correspond to binding energies of 
the minimum energy hep structure and vertical orange (light) 
bars correspond to the binding energies of the minimum en- 
ergy icosahedral structure. The inset shows the trend in the 
bulk binding energy per atom for Mn, Fe, Co, Ni and Cu. 
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FIG. 6: (Color online) Plot of binding energies of the mini- 
mum energy hep and icosahedral structures for each Fen and 
Co„ cluster (15<n<20) in spin-polarized calculations. The 
squares correspond to the data points for the Fe„ clusters 
(solid squares for the hep symmetry and empty squares for 
the icosahedral symmetry) and the circles for the data points 
of Co„ clusters (solid circles for the optimal hep structure and 
empty circles for the optimal icosahedral symmetry). 



eratc for both spin polarized as well as non-spin polarized 
calculations. 

For Mn clusters, the effect of noncollinearity has been 
discussed in literature Mn is prone to noncollinearity 
due to the presence of competing nature of magnetic in- 
teractions, though the degree of noncollinearity is found 
to decrease for cluster sizes larger than 13.— For Fe and 
Co clusters, the degree of noncollinearity is reported to 
be further smaller compared to Mn.— Noncollinearity is 
favored by the magnetic energy associated with larger 
magnetic moments, which competes with chemical bond- 
ing energy. One would therefore expect reduction of non- 
collinearity in moving to larger cluster size as well as mov- 
ing from Mn to Fe and Co. However, in order to check 
the influence of the possible noncollinearity which may 
arise due to competing magnetic interactions as well as 
orbital component of magnetic moment, driven by spin- 
orbit (SO) coupling, we have repeated our calculations 
for Fei9 and Coig clusters in terms of GGA-fSO cal- 
culations. The obtained results indicate that the Feig 
and Coig clusters are essentially coUinear, with degree 
of noncollinearity being less than 2° in agreement with 
that reported previouslyi^ Though the orbital compo- 
nents of magnetism are found to be of finite values « 0.08 
/is, importantly the calculations carried out considering 
noncollinearity leads to only small changes in the bind- 
ing energy differences of the icosahedral and hexagonal 
geometries by 1- 2 %, keeping the main conclusion of our 
study unchanged. In the following, we focus primarily 
on the Fcig and Coig clusters, for which the switching of 
the stable phase between hep and icosahedral structures 
occurs. 

It is important to note that the trend in binding en- 



ergy calculation is very robust, being independent of the 
type of pseudo-potential or the nature of the exchange- 
correlation functional used in this study. We also found 
that this trend to be general for clusters having sizes 
15<n<20^ The structures for n=15, 16, 17, 18 and 20 
were obtained by removing or adding atoms from the op- 
timal 19-atom cluster structure and then letting them to 
optimize for all possible collinear spin configurations of 
the constituting atoms. In Fig. [6j we have shown a plot 
of binding energies of the optimal hep and the optimal 
icosahedral structures of the Fe„ and Co„ clusters consid- 
ering the clusters sizes in the range 15<n<20. It clearly 
indicates that the icosahedral growth pattern is more fa- 
vorable for the small Fe„ clusters and the hep growth 
pattern for the small Co„ clusters in the spin polarized 
calculations, in agreement with the trend observed for 19 
atoms clusters discussed above. In Table II, we have also 
shown our estimated magnetic moments of the optimized 
hep and icosahedral structures of Fe„ and Co„ clusters 
in this size ranges. Notice that our estimated magnetic 
moments for the optimized structures are in fair agree- 
ment with the recent result of Stern-Gerlach experiments 
for Fe-clusters^ and Co-clustersi^ 

In order to understand the optimal structures and the 
distortions in the structure that arises during optimiza- 
tion procedure, in Table III, we list the rms distortion of 
the bondlengths in the optimized geometries, which gives 
us the feel of the distortions that accompany optimiza- 
tion. 

The pertinent question, therefore, is what drives this 
phenomena ? To see the effect of magnetism, we first 
calculated the magnetic energy which is defined as the 



TABLE II: Calculated magnetic moments of the optimized 
hep and icosahedral structures of Fe,i and Co„ clusters 
(15<n<20) in spin-polarized calculations. For comparison, 
we have also listed the recent experimental values (Ref. 59 
for Fe„ clusters and Ref. 60 for Co„ clusters) of magnetic 
moments in this size ranges. 



Clusters 


M {fiB 


/ atom) 




Clusters 


M [flB 


/atom) 




Theory 


Expt. 




Theory 


Expt. 




hep 


leosa 






hep icosa 




Fei5 


3.07 


3.20 


2.72 


Coi5 


2.07 1.93 


2.38 


Fei6 


3.13 


3.13 


2.94 


Coi6 


2.13 1.88 


2.53 


Fei7 


3.18 


3.06 


2.86 


Coi7 


2.06 2.06 


2.24 


Feis 


3.11 


3.11 


3.02 


Coi8 


2.00 2.00 


2.07 


Fei9 


3.05 


3.05 


2.92 


Coi9 


2.05 1.95 


2.21 


Fe2o 


3.00 


3.00 


2.73 


C020 


2.00 1.90 


2.04 



0.8 
I 0.6 

tsf 0.2 





■ ■ hep 


(b) . 


0.7 




icosa 




0.6 






\ ; 






"^^^■« . 


0.5 










1 1 1 




1 1 - 


0.4 



Mn^j, Fcjg Cojg Nii9 Cuj^ Fe^g Co^^ 



FIG. 7: (Color online) (a) Magnetic energy (Em), calcu- 
lated as the energy difference between the spin polarized and 
non-spin polarized calculations of the optimized hep (solid 
squares) and icosahedral (light solid dots) structures plotted 
for each X19 cluster (X — Mn, Fe, Co, Ni, Cu). (b) Zoomed 
plot around the data points for the Fei9 and C019 clusters. 



TABLE III: Calculated root mean square (rms) deviations of 
bond lengths of Co,i and Fe„ clusters for n= 15 - 20, for both 
hep and icosahedral symmetry based optimized structures. 
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energy difference between the magnetic (spin-polarized) 
and nonmagnetic (non-spin-polarizcd) calculations for 
each of the hep and icosahedral structures of X19 clus- 
ters, estimated for their optimal structures in magnetic 
and nonmagnetic calculations. Fig. [7{^a) shows the plot 
of magnetic energies of the X19 clusters for the hep and 
icosahedral symmetry based optimal structures. It is in- 
teresting to note that the magnetic energy of icosahedral 
structure is much higher than that of the hep symmetry 
based structure for the Feig clusters (and also for Mnig 
cluster, though we do not bring it into our discussion 
due to the assumption of collinearity in our calculation 
as mentioned before). On the other hand, it is of similar 
magnitudes for the hep symmetry based structure of C019 
cluster and its icosahedral counterpart, with hep being 
somewhat higher. The magnetic energy difference be- 
tween hep and icosahedral structures is negligibly small 
in case of the Niig and Cuig clusters. We have shown 
the zoomed plot again around the Feig and C019 data 
points in Fig. ^h) in order to see the effect of magnetic 
energy more closely, which shows opposite trend of mag- 
netic energy gain between Feig and that of Coig more 
clearly. We note that the difference of magnetic energy 
gains between the hep and icosahedral structures in the 
case of the Coig cluster, is relatively small compared to 



that of the Fcig cluster. As the d-shell gets progressively 
filled up, starting from the half-filled situation with the 
highest atom-centered magnetic moment for Mn„ clus- 
ter, the magnetic energy gain gets progressively weaker, 
thereby the role of magnetism being more important for 
Fe, compared to Co. 

In order to understand the gain in magnetic energy for 
the icosahedral structure of the Fcig cluster and for the 
hep structure of the Coig cluster, we have studied the 
density of states (DOS) of the optimized hep and icosa- 
hedral structures of the Feig and Coig clusters for both 
the magnetic and nonmagnetic calculations as shown in 
Fig. [51 We note that compared to the nonmagnetic DOS, 
the gap in the majority spin channel is significantly en- 
hanced in case of the icosahedral structure of Fcig and 
the hep structure of Coig, indicating their enhanced sta- 
bility. On the other hand, for the optimal hep structure 
of Feig and for the optimal icosahedral structure of Coig 
cluster in case of spin-polarized calculation, there are fi- 
nite amount of states around the Fermi energy, which 
reduces the stability of the system, compared to that of 
the corresponding icosahedral and hep structures. 

We next study another relevant quantity which has 
been used previously to examine the relative stability 
between the various classes of isomers for the 3d late 
transition metal clusters, namely the hybridization of the 
atomic 3c? and 4s orbitals. The s-d hybridization index as 
quantified by Hakkinen et aS^ and later used by Chang et 
as well as Wang et for transition metal clusters, 
is defined for a 19 atoms cluster as 
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where w( ^ {w( ^) is the projection of i-th Khon-Sham or- 
bital onto the s (d) spherical harmonic centered at atom 
/, integrated over a sphere of specified radius. The spin 
index is implicit in the summation. Our calculated s-d 
hybridization index for the optimized structures of both 
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FIG. 8: (Color online) Plot of DOS of the optimized hep and icosahedral symmetry based structures of Feig cluster (upper 
panels) and C019 cluster (lower panels) in spin-polarized as well as non-spin polarized calculations. The results for nonmagnetic 
calculations have been shown with shaded curves. The smearing width is fixed at 0.1 eV. The vertical line through zero is the 
Fermi energy for each system. The inset shows the DOS corresponding to the spin-polarized calculations, zoomed around the 
Fermi energy with a smearing of 0.001 eV. 



the symmetries for the Fcig, C019 and also for the Niig, 
Cuig clusters have been plotted in Fig. |9l To see the 
effect of magnetism, we have studied the s-d hybridiza- 
tion of the optimized structure of each cluster for both 
the magnetic and nonmagnetic phases. 

It is seen that the optimized hep structures have con- 
sistently higher values of Hsd than those of the optimized 
icosahedral structures of 3d late transition metal clusters 
Xig for both spin-polarized and non-spin polarized cal- 
culations. In order to see distinctly the effect of magne- 
tization on s-d hybridization, we have plotted the differ- 
ence of s-d hybridization indexes between the optimized 
hep and the optimized icosahedral structures for both 
the magnetic and nonmagnetic calculations in the right 
panel of Fig. |9l The positive (negative) value of this 
difference, Ksd'^^^^ , indicates that the hep (icosahedron) 
structure has higher s-d hybridization. It is clearly seen 
that though the difference is positive for all the late tran- 



sition metal clusters, it shows some variation across the 
series. The s-d hybridization gain in favor of hexagonal 
structure is the maximum for the magnetic Coig clus- 
ter, showing a factor of about 6 times enhancement com- 
pared to nonmagnetic Coig. Cuig cluster being essen- 
tially nonmagnetic, s-d hybridization gain between the 
two structural symmetries remains same both in mag- 
netic and nonmagnetic calculation of Cuig. The s-d hy- 
bridization gain remains similar for the magnetic Feig 
and nonmagnetic Feig (Hgd'^''^^ ~0.3) and that for mag- 
netic Niig and nonmagnetic Niig {Hgd'''^^^ '^0.4). We 
therefore, conclude that the gain in s-d hybridization sta- 
bilizes the hep symmetry based structure over the icosa- 
hedral symmetry based structure for the Coig cluster. 
This is also helped in a way by the small but positive 
magnetic energy gain in favor of hep phase of the Coig 
cluster. So the s-d hybridization helped by magnetic en- 
ergy gain stabilizes the hep symmetry based structure in 
case of the Coig cluster. On the other hand, for Feig 
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cluster, the large magnetic energy gain in favor of the 
icosahedral symmetry, decides the final stability, thereby 
counteracting the hybridization energy gain in favor of 
hexagonal symmetry. 



IV. SUMMERY AND CONCLUSIONS 

To summarize, we have investigated the relative stabil- 
ity of the 3d late transition metal clusters specially of 19 
atoms between hep and icosahedral symmetries. Among 
all the members, the Coig cluster prefers an unusual sta- 
bilization in hexagonal symmetry, while the rest show 
the preference of icosahedral symmetry. Our study nicely 
demonstrates that this curious result is driven by the in- 
terplay of the gain in magnetic energy vis a vis the gain 



in hybridization energy. For the Coig clusters, the hy- 
bridization energy gain helped by magnetic energy gain 
favors the stabilization of hexagonal symmetry, while for 
clusters like Feig, the large magnetic energy gain in icosa- 
hedral symmetry topples the s-d hybridization gain in fa- 
vor of hexagonal symmetry and stabilizes the icosahedral 
phase. We find the obtained trend to hold good as well 
for clusters having sizes between 15 and 20. 
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FIG. 9: (Color online) Plot of (a) s-d hybridization index for the optimized hep and optimized icosahedral 19-atoms clusters 
of Fe, Co, Ni and Cu, both in the spin-polarized (magnetic) calculation and non-spin polarized (nonmagnetic) calculation 
(solid circles and squares represent the results for spin-polarized calculation and empty circles and squares correspond to the 
data points for non-spin polarized calculations), (b) the difference in s-d hybridization (Rsd'^^^ ^ =Ksd'^'^^-'ilad"^°"'^) between the 
optimized hep and icosahedral structures of 19-atom clusters both in the spin-polarized and non-spin polarized calculations. 



